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ABSTRACT 

We generalize the predictions for the CMB anisotropy patterns arising in Bianchi type VI L, universes 
to inclu de a dark energy component. We consider these models in light of the result of l.Taffe et al. I 
d2005albl) in which a correlation was found on large angular scales between the WMAP data and the 
anisotropy structure in a low density Bianchi universe. We find that by including a term Ha > 0, the 
same best-fit anisotropy pattern is reproduced by several combinations of cosmological parameters. 
This sub-set of models can then be further constrained by current observations that limit the values 
of various cosmological parameters. In particular, we consider the so-called geometric degeneracy in 
these parameters imposed by the peak structure of the WMAP data itself. Apparently, despite the 
additional freedom allowed by the dark energy component, the modified Bianchi models are ruled out 
at high significance. 

Subject headings: cosmology: cosmic microwave background - cosmology: observations 


1. INTRODUCTION 

While cosmology appears to be converging on a “con¬ 
cordance model” that describes the universe as infla¬ 
tionary and isotropic, there remain unexplained anoma¬ 
lies in the CMB data, and other models are not yet 
ruled out. The WMAP data provide some of the most 
accurate measurements yet of the cosmic microwave 
background and contribute to hig h accu racy de termi¬ 
nations of cosm ological parameters iBennett et a.ll[2003t 
iSnergel et al. 112003 1. However, there are several stud¬ 
ies that show that at large scales, the CMB is not sta¬ 
tistically isotro pic and Gaussian, as pred icted by infla¬ 
tion ■ theor^_jd^Oliveir^Cost^t^2|20(2[Eriksei^OiD 
120041 iHanse^^t ^hll200^1WeivaetaLll20Mb 

In [Jaffe et al. (l2005alH) . we examined a particular set 
of anisotropic cosmological models: the Bianchi type 
VII/, class of spatially homogeneous generalizations of 
Friedmann universes that include small vorticity (uni¬ 
versal rotation) and shear (differential expansion) com¬ 
ponents. Surprisingly, we found evidence that one of 
these models correlates with the CMB sky and that sub¬ 
tracting this component resolves several of these observed 
an omalies. The models used in that study were derived 
bv lBarrow et al. I ( 1985 ) and include no cosmological con¬ 
stant component in the total energy density. The best-fit 
model found by that study required fi m o = 0.5, implying, 
H t. = 0.5, i.e., a sign i ficant ly negatively curved universe. 

lLand fc Magueiiol ( 2005 ) subsequently considered flat 
Bianchi models and sought explicitly to resolve the prob¬ 
lems of the low quadrupole and the low-t? alignments us¬ 
ing a statistic constructed to achieve that purpose. Their 
analysis does indeed find a model that fixes these anoma¬ 
lies, but it remains statistically insignificant as a detec¬ 
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tion. Our original result has the benefit of being a de¬ 
tection that is based entirely on a simple least-squares 
fit of the Bianchi template to the data, yet we find that 
it serendipitously resolves several anomalies without that 
requirement having been built-in to the search algorithm. 
The importance of our result, therefore, lies in the fact 
that it resolves the problems of the low quadrupole, low- 
11 alignments, power asymmetry, and non-Gaussian cold 
spot. 

However, both analyses neglect the fact that the exist¬ 
ing Bianchi solutions include no dark energy component. 
Furthermore, the best-fit results require an energy con¬ 
tent that is consistent neither with other astronomical 
observations nor with the CMB itself. In this work, we 
present a modification of the Barrow et al. Bianchi type 
VII/, solution so that it includes a cosmological constant 
term in the evolution. We discuss the impact of that 
term on the structure of the resulting CMB anisotropy 
pattern, particularly the degeneracy that is introduced 
in the model space by the addition of Ha- We then dis¬ 
cuss the viability of A models that are morphologically 
identical to our best-fit Ha = 0 models by considering 
the constraints imposed by different measurements of the 
cosmological parameters. 

2. BIANCHI MODELS WITH A > 0 
2.1. Solution 

Firs t we shall generalize the equations in lBarrow et al. 1 
(Tl985l to include a cosmological constant A. The ba¬ 
sic assumption is that the type VII/, universe is close to 
FRW. 

We start with the equations of motio n usin g expansion- 
norm alized variables (see, e.g., IWainwright V Elli4 
4997 ). The equations describing the evolution of Bianchi 
type VII/, u niverses with a tilting perf ect fluid can 
be fou nd in IColev V Hervikl ((2005) and iHervik et aiJ 
(l2005al) . In our case, we consider a tilting perfect fluid 
with 7 = 1 (dust) and a co s molog ical constant. The 
equations in IColev fc Hervikl ( 2005 ) can easily be gen¬ 
eralized to include a cosmological constant by adding a 
matter term (Ha, 7a = 0) wherever there is a matter 
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term (f 2 , 7 ). 

Furthermore, we assume that the universe is close to 
a FRW at all times and that the tilt velocity of the dust 
fluid is small; hence we assume: 

E 2 <1, |N x | <C 1, U!<1, |v| < 1. (1) 

Here, E 2 = a 2 /(3H 2 ) is the expansion-normalized shear 
scalar, N x is a complex curvature variable, and rq, v = 
V 2 +W 3 are the tilt components (see lColev fc Hervikl2005l 
for details). The curvature variables A and N need not 
be small, and they are related to the group parameter h. 
as follows: 

A 2 ~ 3hN 2 . 

We also assume that the parameter h is not too small: 
h > 0(1). (As h —> 0, the assumptions made for this 
derivation break down and a qua litatively different so lu- 
tion for type VIIo is needed. See lHervik et alJfeOOfiU ') 
Given the above assumptions, the deceleration param¬ 
eter is 


1 . 


:H m Ha, 


( 2 ) 


and the equation of motion for the Hubble scalar is 

H' = -(q+l)H, (3) 

where the prime indicates the derivative with respect to 
the dimensionless time T defined by dt/dT = 1/H 1 and 
t is the cosmological time. 

At the lowest order, the tilt components that are re¬ 
lated to the vortici ty are the component s V 2 + A 3 = v 
(in the notation of iColev fc HerviUl2005l) . These induce 
non-zero shear components E 12 + iSi 3 = Si via the lin¬ 
ear constraint equation, which in the FRW limit reduces 
to: 


Y^iN(i — 3Vh) + fi m v = 0 . 


( 4 ) 


Moreover, close to FRW the curvature K = A 2 and mat¬ 
ter densities fi m and Ha are related via the Friedmann 
equation 


1 — Ha + A + fi m . 
We also define x to be 


An 


1 — Hao — H 


mO 


( 5 ) 


( 6 ) 


Eq. QJ now reduces to Barrow et al.’s eqs.(4.6) and 
(4.7) by dropping the subscript 0 and replacing H with 
H m . The equations of motion can now be solved. 

By choosing T 0 = 0, x = e a °Hn, we can relate z to T 
and the time variable r in iBarrow et al. to the 

redsliift z: 


1 + z = e 
dT 


dr 


= He T e a °. 


( 7 ) 

( 8 ) 


The Hubble scalar, K , fi m , and Ha can then be written 


H(z) = H a H(z), 


N(z) — [Hao + Aq(1 + z)" -(- H m o(l + z) 3 ] 


31 2 


( 9 ) 


( 10 ) 


K = 


H m — 


Ao(l + z) 2 
n 2 (z) ’ 

H m o(l + ~) 3 


Ha = 


n 2 (z) 

Hao 


( 11 ) 


( 12 ) 


(13) 


n 2 (zy 

where HAo,ATo,H m o are the present values of the 
expansion-normalized A, curvature and ordinary matter 
(dust) and obey the constraint: 

1 = Hao + Ao + H m0 . 

We also note that K = fi*, > 0 implies negative cur¬ 
vature. The shear components E 12 and E 13 , and the 
remaining curvature variable are given by 

El2,o(l + ~) 3 


Ei o — 


E 13 — 


N = 


H(z) 

£13,0(1 + z ) 3 

H(z) 1 

N 0 (l + z ) 


H(z) 


Conformal time becomes 


1 [ z dz 

w>' 


(14) 


(15) 


(16) 


(17) 


which can be used in the Barrow expressions for s and 
ip. The integrals of Barrow et al. equations (4.12) and 
(4.13) become 


/•ze 


s(l — s 2 ) sin^(l + z) 2 dz 


Jo (1 + s 2 ) 2 H(z) 

r ZE s(l — s 2 ) cosi/>(l + z) 2 dz 

/ 0 


(18) 


(19) 


(1 + s 2 ) 2 H(z) 

with C 3 = 4, where ze is the redshift at photon emission. 
In calculating the constant Ci, the density Ho should now 
be fi m o, i.e. the matter density only. 

The following assumptions should be kept in mind re¬ 
garding this solution: 

• The tilt velocity is small. In a dust- or A-dominated 
Bianchi type VII^, universe, the tilt asymptotically 
tends to zero at late times. Using the above re¬ 
lations, it can also be shown that for the case of 
interest, it remains small back to the last scatter¬ 
ing surface as well. 

• The universe is close to FRW at all times. 
This assumption is required to make the problem 
tractable. In a A-dominated universe, the type 
VIU models become close to FRW. For the cases 
in question, plugging the numbers into the above 
relations shows that this assumption still holds at 
the surface of last scattering. (At very early times, 
however, this is not the case.) It is likely that mod¬ 
els for which this assumption does not hold at last 
scattering would then have greater shear effects, 
but the details would require a more complicated 
derivation. 
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• iBarrow et al. I ll198511 only consider the contribu¬ 
tion from the vorticity that directly involves the 
shear variables £12 and £ 13 . The (expansion- 
normalized) shear scalar, £ 2 = a 2 /(3H 2 ), is given 
by 

£ 2 = £^ + |£ x | 2 + £? 2 + £ 2 3 ( 20 ) 

an d involves ther e fore a ll of the shear components. 
In IBarrow et al. I (1985) and here, the other shear 
components have simply been ignored in the cal¬ 
culation of A T/T. Considering additional shear 
degrees of freedom is in principle possible but in 
practice difficult. The effect of such additions is 
one of the largest unknowns in this analysis. 

2 . 2 . Properties 

In the presence of the cosmological constant, a universe 
that is very nearly flat now and that was very nearly flat 
at the time of last scattering can become significantly 
negatively curved at intermediate redshifts. We can see 
this from Eqtn. im above. 

Figure U shows how K(z) evolves for different values 
of S 2 ao and fl m o- At high redshift, the matter dominates 
and the curvature vanishes, and at low redshift, can 
dominate, depending on its exact value. At intermediate 
redshifts, however, the universe may become more nega¬ 
tively curved. For a currently almost flat universe with 
parameters close to the observed values of Dao ~ 0.7 and 
fl m o ~ 0.3, the curvature is at most only a few percent 
(black curve). Only for very small matter densities and 
very large dark energy densities does the curvature be¬ 
come large at intermediate redshifts. To reproduce the 
observed asymmetry with these A > 0 models requires 
either a large current negative curvature (pink and blue 
dot-dashed curves) or a very small matter density (green 
dashed curve). 



Fig. 1.— Examples of the evolution of K(z). See text. 

The addition of a A component adds a third parameter 
to the Bianchi parameter space, but there is a degener¬ 
acy to the resulting morphology of the Bianchi-induced 
pattern. The same structure in the induced anisotropy 
pattern can be reproduced by different parameter combi¬ 
nations. The addition of dark energy results in a tighter 


spiral for a given value of x. because the redshift to the 
surface of last scattering corresponds to a larger time 
difference. The geodesics have therefore completed more 
rotations since recombination for a given value of x when 
Ha > 0. Secondly, the focusing of the spiral depends on 
the curvature between the observer and the last scatter¬ 
ing surface. Significant focusing can arise from a universe 
that is very open now, or from a universe where A is 
large enough to give negative curvature at intermediate 
redshifts (as shown in Fig. [T]) 

Three combinations of parameters, for example, can 
reproduce the same apparent structure on the sky : 

1. (a;, S2 ao, ^mo) = (0.62,0, 0.5) - i.e., no A and a mat¬ 
ter content half critical, giving a large cur rent neg¬ 
ative c urvature (the model found in I.Taffe et al. I 
I2005H1 : 

2. (a;, flAOj ^mo) = (0.8, 0.4,0.3) - i.e., approximately 
the observed matter content (baryon plus dark) , 
with some A but still a large current negative cur¬ 
vature; 

3. (x, Dao, ^mo) = (4.0,0.987,0.003) - i.e., a much 
smaller than observed matter content, but a nearly 
flat current curvature (S\ = 0.01); the large A 
causes a large negative curvature at intermediate 
redshifts. 

The real difference between these models is the ampli¬ 
tude: has an amplitude ~ 80% of 0 , and 0 has 

an amplitude of ~ 8 times that of (1). The shear and 
vorticity values obtained from fitting these templates to 
the data would change accordingly. 

Effectively, this implies that a given structure on the 
sky characterized by the amount of geodesic focusing and 
the number of spiral turns is the same for all models 
along a line in the three dimensional parameter space of 
(fiko.^AO,*)- The degeneracy in the templates is only 
broken by the amplitude of the variation. But (a/H) 0 is 
what we measure by fitting a template to the sky, so we 
cannot distinguish among the degenerate models without 
an independent measurement of the shear. 

It is not straightforward to calculate where such lines 
of degeneracy lie. Instead, we determine this empirically 
by generating a grid of models and comparing them to 
the previously determined best-fit model. 

3. VIABILITY OF A BIANCHI TYPE VII w MODEL 

We now address the question of whether the degener¬ 
acy of the Bianchi type VII/i models with Hao > 0 de¬ 
scribed in 0gives us enough freedom to define a Bianchi 
model th at is mor pholog ically identical to the best-fit 
model in I.Taffe et al. I d2005aiq ) and that is also consis¬ 
tent with estimates of the cosmological parameters. 

3.1. Matter Density and ’Geometric’ Degeneracy 

lEfstathiou fe Bondi dl999ft describe the limitations of 
constraining cosmological parameters with CMB data 
alone. In particular, there is a ’geometric’ degeneracy in 
the curvature and dark energy density (or equivalently, 
the matter density.) Degenerate models have: 1) the 
same values of the physical baryon density , Wb = h, 2 , 

and cold dark matter density, u c = 2 ) the same 
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Fig. 2.— Degeneracy contours on the Da — D*. plane for Ct>b = 
0.023 and a) m = 0.13. Along the red line denoting TZ for the WMAP 
best-fit parameters, the circles mark example power spectra that 
are shown in Figure 0 


primordial fluctuation spectrum; and 3) the same value 
of the parameter 


n = 


{ sinh [y/uJky\ if u k > 0, 

y/UkU if Wfe = 0, 

sin[y AJ^y\ if to k < 0, 


where 


da 


V = 


y/bJ m a + IVkO 2 + WA« 4 


( 21 ) 


( 22 ) 


(where a rec is the scale factor at recombination). These 
models will have almost the same power spectra. Con¬ 
tours of constant TZ in model space flk versus are 
shown in Figure |21 with red highlighting the curve 
intersecti ng th e best -fit, flat, WMAP-only parameters 
(iSnergel et al. 11200 11 Table 1). Figure 0 shows power 
spectra for those degenerate models (calculated using 
CMBFAST. 4 ) When normalized, the spectra are iden¬ 
tical at the acoustic peaks, and the degeneracy only fails 
at large angular scales. 

In the three dimensional parameter space of 
(flfcOi f^AO) x), the geometric degeneracy forms a sur¬ 
face, while the Bianchi models with identical structure 
form a line. The two need not intersect, but in the 
case of the observed WMAP CMB-only power spectrum 



Fig. 3.— Power spectra for several models degenerate with the 
WMAP data (marked with circles in figure [2]) The solid versus 
dashed green and blue curves are where there are two degenerate 
models at the same value of the dashed is that with the higher 
The spectra are normalized so that the amplitude at the first 
acoustic peak is always the same. The dotted black line shows the 
cosmic variance uncertainty for each data point around the WMAP 
model. 


of iSoergel et al. I (1200311 and best-fit Bianchi model of 
l.Taffe et al. I (l2005alblf . they do. There is a Bianchi model 
that has the identical structure to the best-fit I2 a = 0 
model and lies on the WMAP degeneracy curve. It has 
parameters (flfco,^ ao, x) = (0.028,0.96,2.5). But how 
viable is this region of parameter space? 

As shown in Figure^ the geometric degeneracy is bro¬ 
ken at large angular scales. Models with high ST\o and 
low fl m0 predict too much large scale power. WMAP 
data alone place relatively loose constraints on S2 m o, but 
even these rule out such low values as required for the 
Bianchi models. (The power spectrum of the sky cor¬ 
rected for the Bianch i componen t has less large-scale 
power overall (l.Taffe et al. I2005alblf than the uncorrected 
power shown in that figure, so the problem will only be¬ 
come worse.) 

To quantify these limits, we examine the posterior 
likelihood from t he WMAP data alone using the COS- 
MOMC 5 code of Ibewis Rr. Bridlel 1 2003 1, which imple¬ 
ments a Monte Carlo Markov Chain method. For a 
simple look to compare the WMAP constraints to the 
Bianchi degeneracy, we explore the five-parameter space 
of: 9 (the ratio of the sound horizon to the angular diame¬ 
ter distance), r (optical depth), f Ik (curvature), n s (spec¬ 
tral index of scalar perturbations), and A s (amplitude of 
scalar perturbations), with fixed parameters u}\, = 0.023 
(the physical baryon density), Wdm = 0.107 (the physical 
dark matter density), /„ = 0 (neutrino fraction of dark 
matter density), w = — 1 (dark energy equation of state), 
r = A t /A s = 0 (ratio of tensor to scalar fluctuations). 
Other parameters such as Hq, and f2 m o are derived 
from this set. (The free parameters have only such pri¬ 
ors as defined by ranges that are much broader than any 
realistic uncertainty.) The resulting constraints in the 
— Ha plane (marginalizing over the other parame¬ 
ters) are shown in Figure El in red. All Bianchi models 
with the structure of the best-fit model lie well outside 
the 95% confidence region from the WMAP data alone. 


4 http://www.cmbfast.org/ 


5 http://cosmologist.info/cosmomc/ 
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0 

Fig. 4.— Par ameter space by f2 m . In green are the super¬ 
novae constraints JKnop et al. 1120031) . The solid red line is the 
same geometric degener acy curve shown in Fig ure l2l In blue are 
the HST Key Project IFreedman et aid 1200U) constraints from 
O m o = ^m/h 2 (assuming fixed o; m and u c )- The black solid lines 
show an approximate representation of where the Bianchi degen¬ 
eracy lies, from the original model at (x, ^mo) = (0.62, 0, 0.5) 
to the one that lies on the WMAP geometric degeneracy curve at 
(0.028,0.96,2.5) (each shown with a violet diamond.) The likeli¬ 
hood contours from WMAP data alone (computed using cosmomc; 
see text) are shown with the red dashed lines. In solid magenta 
are the contours from WMAP, supernovae, HST, and SDSS data 
combined, where u; m and u c also vary. 



1 10 100 1000 
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Fig. 5.— Power spectra as a function of optical depth, r, 
each normalized such that the first peak heights are all the same 
(effectively changing A s .) The top panel is the TT power, and 
the bottom the TE cross-power. Lowering r lowers the large scale 
power, but not enough to bring the high Q ,a models down to the 
range of the data. 


Allowing the physical matter densities to vary as well but 
including other datasets, the constraints (in magenta) are 
e ven tighter. 

ITegmark et al. (12004H give D m0 = 0.57+° ^ from 
WMAP data alone, which tightens to D m0 = 0.40^q’o 9 
adding da ta from the Sl oan Digital Sky Survey (SDSS). 
See also iSanchez et al. I ( 2005 1 for CMB+2dFGRS re¬ 
sults. 

3.2. Optical Depth 

Figure 0 shows that the geometric degeneracy breaks 
at large angular scales, and the models with high Da 
have too much large scale power. But there is an ad¬ 
ditional degeneracy if we allow the optical depth, r, to 
vary. Adjusting r and A s such that A p = A,e _2r remains 
constant has exactly the effect we need of modifying only 
the large scale power while lea ving the peak heights un¬ 
changed ITegmark et al. l[2004il . Figure 0 shows how this 
works, but that the effect is not large enough. Even a r 
of zero does not bring the power at large angular scales 
within range of the data, and that model is inconsistent 
with the large-scale peak in the TE spectrum. 

One can also ask if the addition of the Bianchi compo¬ 
nent can affect the variance of the large scale TE cross¬ 
power, resulting in an incorrect estimate of r due to 
a chance alignment of the polarization signal with the 
Bianchi structure. Simulations with no reionization and 
with an added Bianchi temperature pattern (the same lo¬ 
cation and shear amplitude as our best-fit model) show 
that this is not the case. The cross-power remains flat at 
low £ with t he expected va r iance. 

Note that lHansen et al. d2004a l find different values 
of r derived from fitting temperature data in the North¬ 


ern and Southern hemispheres as defined in the frame 
of reference t h at ma ximizes the power asymmetry (see 
lHansen et a,l.ll2004hl. l The additional large-scale struc¬ 
ture in the South that the Bianchi template reproduces 
could cause this effect. If a polarization signal were also 
produced in a Bianchi geometry and correlated with the 
temperature anisotropy, it would influence the measured 
optical depth from the low-£ TE peak as well. 

3.3. Supernovae, H 0 , etc. 

The degeneracy contours in Figure 0 are plotted 
for co nstant values of = 0.13 and ccb = 0.023 
(from iSoergel et al. 11200311 . For different points in the 
(Dfco, Dao) space, D m o changes and therefore the value of 
h = Hq/(100 km s _1 Mpc _1 ) changes, increasing with in¬ 
creasing Dao along the degeneracy curve. For models of 
high Dao that lie on the geometric degeneracy curve, Hq 
reaches values over 300. This is rule d out at h i gh signif¬ 
icance by the WMAP data itself; ITegmark et al. u21)l) 1: 
give Hq = 481^2^111 s _1 Mpc _1 . Independent determi¬ 
nations of Hubble’s consta nt also rule out these v alues, 
e.g., the HST Key Project llFreedman et al. ll2001fl value 
of 72 ± 8km s -1 Mpc -1 . See, e.g., iSnergel et al. I <120031 
Table 3) for other estimates of the Hubble constant. 

High-redshift supernovae observations can break the 
geometric degeneracy by independently placing con¬ 
straints that are nearly perpendicular to the WMAP 
constraints in D m — Da space (see Figure 0) The de¬ 
generacy curve for the best-fit Bianchi model compared 
to the WMAP data lies outside the 98% confidence con¬ 
tour determ ined by the Supernova Cosmology Project 
( Knoo et al . 120031. fit no. 6). 

It is worth pointing out that it is not clear whether 
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type la supernovae are truly standard candles out to high 
redshift. The uncertainty, however, is not enough to ac¬ 
commodate such low values of O m o- 

The addition of cosmic vorticity and shear would, of 
course, influence the determinations of such parameters. 
The sky coverage of supernova, for example, is fairly ex¬ 
tensive outside the Galactic Pla ne region, but the sam- 
ple size is small. Studies such as lKolatt & Lahav 1200 if) 
have ruled out significant dipole or quadrupole asymme¬ 
try in the expansion, but a test for a more complicated 
anisotropy structure induced by the vorticity would re¬ 
quire far more supernova in the sample and a good sky 
coverage. Furthermore, as the current value of the shear 
expansion is very small, distance measurements to rela¬ 
tively low red-shift (compared to the CMB) objects may 
not be sensitive enough to detect it. This also implies 
that such a small current shear would not significantly al¬ 
ter measurements of other cosmological parameters such 
as the Hubble constant that are dependent on relatively 
low-redshift observations. 

3.4. Other Bianchi Models 

One can also ask, then, how much one can vary the 
parameters of the Bianchi model and still have a statisti¬ 
cally significant fit. The most interesting region on Fig- 
ureQ]is near (O m ,OA) = (0.15,0.75), where the WMAP 
and supernovae constraints approach the region of the 
best-fit Bianchi model. This model resembles the best- 
fit model, but since it is closer to flat than the Bianchi 
degeneracy curve, there is less geodesic focusing. As a 
result, the model is not so good a fit to the data, and 
has a significance, compared to Gaussian realizations, of 
only 85%, compared to over 99% for the model at the 
same but on the Bianchi degeneracy curve. 

Considering that the WMAP data somewhat favor a 
closed universe, one might ask about closed models with 
vorticity and she ar. These are Bi anch i type IX models, 
also discussed in iBarrow et al. Il98l . Unlike the open 
type VIU models, however, closed models exhibit neither 
geodesic focusing nor the spiral pattern, even in the pres¬ 
ence of vorticity. Barrow et al. use them to place limits 
on vorticity and shear simply using the quadrupole. Such 
closed models do not reproduce the morphology needed 
to explain the power asymmetry or the cold spot, and 
it would be impossible to claim any detection with only 
the quadrupole as an observable. 

3.5. Other Dark Energy Models 

The cosmological constant is the simplest form of dark 
energy, and all observations so far remain consistent with 
it. Other models are not ruled out, however, and in 
some cases are favored. Alternatives come in many vari¬ 
eties, some with physical motivation, ot hers constructe d 
to gi ve a particular result. (See, e.g., iPadmana bhan 
f2003h Here, we address the question of whether an alter¬ 
native dark energy model can bring our Bianchi pattern 
any closer to the constraints imposed by the data. 

Dark energy models are characterized by their equation 
of state, p = wp, where a cosmological constant A corre¬ 
sponds to a model with a constant w = —1 . Alternative 
theories allow w to vary with time, as in “quintessence” 
and “k-essence” models. A rather ad hoc parameteriza¬ 
tion is often used of the form w = wq + w\Z, which allows 


compariso n of generic dark energy models with s uper- 
nova data <|Wang et alJl2004t iDicns fc Renkoll20f)4l) . 

The derivation in m can be generalized for any dark 
energy model by the appropriate substitution into equa¬ 
tions m and d of fix instead of Da- We then have 

H(z) 2 = Cl X of x (z) + K 0 (l + z ) 2 + n m0 (l + zf (23) 

where the function fx is derived from the equation of 
state as 

fx=ex p| 3 ^ Y^j (i +w(z'))} ( 24 ) 

Examples for extreme values of Wq and W\ are shown 
in Figure El As was shown in Figure ^ we can look 
at the evolution of the curvature as a function of time 
for these models and see if any display the required 
negative curvature at intermediate redshift. Using this 
parameterization, or alternatively the parameterization 
w = Wo + ivi-^p [’ does not give much additional free¬ 
dom, however, to create a viable model. Ultimately, the 
amount of curvature and therefore focusing is still largely 
a function of the relative densities of matter and dark 
energy and of the current curvature. Manipulating the 
evolution of the dark energy equation of state simply 
changes the time scales on which the transitions between 
the different components of the density occur. 



Fig. 6.— Examples of the evolution of K(z) for different dark 
energy models where the current curvature is constrained to be 
small (Ko = 0.01). f2 m o = 0.29 for all except the green triple-dot- 
dashed curve, where fi m o = 0.003. This is the curve shown in Fig. 
IT] with a small matter density, high f2 a , and approximately the 
amount of curvature needed to create the best-fit Bianchi model. 

3.6. Small-scale Structure 

We have so far ignored an additional issue with these 
models related to th e stocha s tic com ponent of the CMB. 
The analysis in l.Iaffe et al. I l|2005albD assumes that the 
observed CMB anisotropies consist of two independent 
signals: the predicted Bianchi pattern, and a stochas¬ 
tic, statistically isotropic component. The latter may be 
generated via inflation or another mechanism, but if it is 
statistically isotropic at the surface of last scattering, it 
might no longer be statistically isotropic when observed 
after traveling through a Bianchi universe. 

The power spectrum at small scales might deviate from 
predictions due to additional structure (dependent on 
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orientation relative to the preferred axis) introduced by 
propagation in a Bianchi metric. We n ote t hat the asym¬ 
metry o bserved in lEriksen et alJ (I2004T) and lHansen et all 
(2004b) extends to £ ~ 40, while the Bianchi model 
has structure only up to t ~ 15. Geodesic focusing 
might cause such a power asymmetry, though it might 
also require that the asymmetry continue to the small¬ 
est ang ular scales, which is not observed. lHansen et alJ 
<l2004bli also find that some unexplained outliers in the 
WMAP power spectrum are associated with different 
hemispheres. 

Furthermore, the generation of the fluctuations at 
the last scattering surface might also be affected by 
the anisotropy, though in the case in question, that 
anisotropy remains very small (< 10 -4 ) at z ~ 1000. 
Essentially, we expect that the Bianchi models could 
also be constrained by the lack of deviations in the 
power spectrum from the best-fit WMAP /concordance 
model, although detailed predictions for such deviations 
are needed. 

The consistency of the acoustic peak scale with the flat 
concordance model would be very difficult for our best-fit 
Bianchi model to explain. 

3.7. Early Universe 

Considering the success of Big Bang Nucleosynthesis 
in explaining light element abundances, any model that 
significantly changes the physical processes at those early 
times must be ruled out. Several studies have thus at¬ 
tempted to place limits on shear expansion and rotation 
by examini ng their effect on the relative Helium abun¬ 
dance, Y. (Barrow 1 ( 1976 4 originally showed that such 
analysis can place strong constraints on shear. The inter¬ 
esting exception, however, was type VIL, models, where 
Barrow found that the CMB remained a stronger con¬ 
straint. Later studies showed that including more com¬ 
plicated e ffects can reverse the tren d of Y with the sh ear 
(see, e.g.. I.Tuszkiewicz et al. lll983|h iBarrow I (11984(1 re¬ 
visited the issue and showed that in some cases, even ex¬ 
tremely anisotropic models may still have the observed 
Helium abundance. Bianchi type VII/j models, ho wever, 
have not been treated in detail since the result of Barrow 
{HU. After demonstrating for several cases, not in- 
cludin g VII/,., that Y increases for anisotropic models, 
iMatzner et al. I ( 1986 1 conjecture that the same trend ap¬ 
plies for all anisotropic models, but they admit that this 
is difficult to prove. 

In short, it is possible that such Bianchi models are 
strongly ruled out by BBN and the observed Helium frac¬ 
tion. But this has not been definitively proved. 

4. DISCUSSION AND CONCLUSIONS 

We have presented solutions for Bianchi type VIK 
type universes that include a dark energy term and ex¬ 
amined how their morphological properties change over 
the parameter space. The addition of dark energy adds 
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a degeneracy such that different combinations of the 
three parameters (f2 m ,flA,a;) can lead to t he same ob¬ 
served st ructure as in the best-fit model of l.laffe et al. I 
(l2005alHh A template can be constructed that has the 
identical structure of that best-fit model and also falls 
on the geometric degeneracy curve for the parameters as 
measured by WMAP data. 

This model, however, lies well outside the 95% confi¬ 
dence region in H m —H a space for WMAP data, ruled out 
by the over-prediction of large scale power. It also lies 
outside the 95% likelihood contours from the Supernova 
Cosmology Project, and is further inconsistent with the 
measurement of the Hubble constant from the Hubble 
Key Project. Bianchi models that are more consistent 
with these other measurements are no longer good fits to 
the WMAP large scale structure. 

One of the most difficult problems for these models is to 
account for the acoustic peak structure. The anisotropy 
at early times might influence the nature of the fluc¬ 
tuations at last scattering, and the geometry could af¬ 
fect the power spectrum on small angular scales due to 
the geodesic focusing between last scattering and the ob¬ 
server. Detailed predictions for these effects would be 
needed, but it is difficult to envision such an anisotropic 
scenario that happened to reproduce the observed acous¬ 
tic peak structure, mimicking the concordance cosmology 
so well. 

There is currently no prediction for the CMB polar¬ 
ization anisotropy in a Bianchi universe, but such a 
geometry-induced signal could provide an additional test 
of these models. If the preferred direction indicated in 
the temperature data is also reflected in the full sky po¬ 
larization data (expected from further WMAP data re¬ 
leases and, eventually, from Planck), there will be even 
more motivation to consider non-standard models. 

We have shown that our best-fit Bianchi type VII/j 
model is not compatible with measured cosmological pa¬ 
rameters, despite the additional freedom from adding 
dark energy. It is worth reiterating, however, that the 
serendipitous discovery of a theoretically derived tem¬ 
plate that correlates well with the data also happens to 
resolve several anomalies that cannot be explained in the 
standard picture. These particular models may not be 
viable, but lacking any plausible scenario for systematics 
or foregrounds to be the source of the anomalies, non¬ 
standard models that reproduce a similar morphology 
merit continued interest. 
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